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Abstract 

The validity of the cosmic no-hair theorem is investigated in the context of Newtonian 
cosmology with a perfect fluid matter model and a positive cosmological constant. It is 
shown that if the initial data for an expanding cosmological model of this type is subjected 
to a small perturbation then the corresponding solution exists globally in the future and 
the perturbation decays in a way which can be described precisely. It is emphasized that 
no linearization of the equations or special symmetry assumptions are needed. The result 
can also be interpreted as a proof of the nonlinear stability of the homogeneous models. 
In order to prove the theorem we write the general solution as the sum of a homogeneous 
background and a perturbation. As a by-product of the analysis it is found that there is an 
invariant sense in which an inhomogeneous model can be regarded as a perturbation of a 
unique homogeneous model. A method is given for associating uniquely to each Newtonian 
cosmological model with compact spatial sections a spatially homogeneous model which 
incorporates its large-scale dynamics. This procedure appears very natural in the Newton- 
Cartan theory which we take as the starting point for Newtonian cosmology. 
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1. Introduction 



Explicit solutions of the Einstein equations can only be found under restrictive assump- 
tions such as high symmetry. The physical significance of these solutions depends on the 
existence of a class of solutions which is general enough to have a direct physical applica- 
tion and whose qualitative behaviour resembles that of the exact solutions. These remarks 
apply in particular to the homogeneous and isotropic models used in cosmology. Thus it is 
of interest to know which properties of these models are stable under small perturbations. 
There is a large literature concerned with linearized perturbations of such models but this 
can at best give a rough indication of the behaviour of small but finite perturbations. (For 
a general discussion of the relations between different concepts of stability see [9].) The 
basic question studied in this paper is the following. Given a homogeneous cosmological 
model, assumed for simplicity to have compact spatial sections, what is the behaviour of 
solutions which develop from initial data which are, in an appropriate sense, close to being 
homogeneous? Do they remain almost homogeneous? Do they perhaps in some sense 
become more homogeneous in the course of time? 

At present no rigorous results exist on the question just raised concerning non- vacuum 
solutions of the Einstein equations with almost homogeneous initial data. It appears that 
the only results on the nonlinear stability of any solution of the Einstein equations are those 
of Friedrich [7] for de Sitter space and Christodoulou and Klainerman [4] for Minkowski 
space. In this paper we study the analogous question for Newtonian cosmological models. 
There are at least two motivations for doing this. Firstly, Newtonian cosmological models 
are of interest for applications since they are used to model structure formation in the 
early universe [15]. Secondly, they may provide useful insights for the relativistic case. 

As matter model we take either an isentropic perfect fluid with polytropic equation of state 
or dust. The models considered expand for ever and have positive cosmological constant A. 
(Some special solutions with A = are also considered for comparison.) The main results 
concern those solutions with isotropic expansion. (The meaning of this terminology will 
be explained below.) In the case of a positive cosmological constant it is shown that the 
solutions corresponding to nearly homogeneous initial data exist for all positive times and 
that the difference between the inhomogeneous and homogeneous solutions tends to zero 
in a strong sense as t — > oo. On the other hand it is shown by an example using dust that 
the quantity [p — p)/p often does not tend to zero. Here p denotes the density and p the 
mean density. Another example using dust shows that when the cosmological constant is 
zero there are solutions which develop singularities in flnite time and which evolve from 
data which are arbitrarily close to homogeneity. 

To put these results in context it is useful to compare with the case of a non-relativistic 
fluid without gravitation. In that case the spatially homogeneous solutions are time inde- 
pendent. It has been shown by Sideris[17] that if the fluid is polytropic rather general small 
perturbations of homogeneous initial data lead to singularities in finite time. Intuitively, 
one of the mechanisms leading to singularities in a fiuid with non-vanishing pressure is 
the formation of shock waves. In the case of dust shell- crossing singularities will often 
occur. The local effects of gravitation cannot be expected to improve this situation and 
may even aggravate it. There is, however, a global effect which can be helpful. This is 
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that the background solution may be expanding and that this expansion may 'pull apart' 
potential shock waves and shell-crossing singularities, thus preventing them from forming. 
Our results show that when A = the expansion is not powerful enough to achieve this but 
when A > it is, at least for small initial perturbations. Another comparison which is of 
some interest is that with Newtonian cosmological models where the matter is coUisionless 
and described by the Vlasov equation. This kind of matter does not have the same kind of 
tendency to form singularities as a fluid and it is possible to prove global existence without 
any restriction on the size of the perturbations[16]. On the other hand less is known about 
the asymptotic properties of the solutions in that case and it is possible that the existing 
estimates for the asymptotic behaviour are not optimal. 

The method of proof of the existence theorem is inspired by the following observation, 
which is of more general interest. Consider the symmetric hyperbolic system 

A^{u)dtu + A\u)diu = (1.1) 

for a vector-valued function u in n space dimensions. One solution of this equation is given 
by w = 0. It is a well known fact (see e.g. [10, 13]) that even arbitrarily small initial data 
for (1.1) often leads to a solution which develops a singularity in finite time. In this sense 
the solution « = is unstable. If, however, this equation is replaced by 

A'^{u)dtu + A'{u)diU + ku = 0, (1.2) 

where A; is a positive constant then the situation changes dramatically. (This is shown ex- 
plicitly for the simple scalar equation ut+uux+u = in [10].) The solutions corresponding 
to any initial data for (1.2) which are small enough exist for all positive time and, in fact, 
converge to zero exponentially as t — > oo. This can be proved using the fact that energy 
estimates for (1.2) show that as long as the solution is small its norm in the Sobolev space 
, s > A, is decreasing and a bootstrap argument. Below a proof of this kind in a situ- 
ation which is slightly more complicated (due to the presence of the Newtonian potential 
term) will be presented in detail. In the cosmological problem the role of the constant k 
is played hy a/a where a is the scale factor of a homogeneous model. For the argument to 
work this quantity must be bounded away from zero (or at least tend to zero sufficiently 
slowly as t — > oo) and it is the positive cosmological constant which ensures this. There is 
also a term representing the perturbed gravitational field which must be controlled. 

The starting point for the study of Newtonian cosmology in this paper is the Newton- 
Cart an theory, which is a geometrical version and at the same time a slight generalization 
of ordinary Newtonian gravitation theory. It is shown that we can split off from the 
full equations a system of ordinary differential equations which describe the evolution of 
the mean density, a fiat metric and a connection on a torus. These are identical to the 
equations which describe the evolution of spatially homogeneous Newtonian cosmologi- 
cal models. This means that it is possible to associate uniquely to each solution of the 
Newton-Cartan theory with compact spatial sections a 'background solution' which is spa- 
tially homogeneous. It is then useful for the analysis of the dynamics to regard the given 
solution as a perturbation of the corresponding background. In the literature on Newto- 
nian cosmology it is usual to split general solutions into the sum of a background and a 
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perturbation. The analysis in this paper clarifies the significance of a splitting of this kind 
and shows that it can be invariantly defined. This contrasts with the situation in general 
relativity where the problem of associating a homogeneous cosmological model with an 
inhomogeneous one so as to give a mathematical formulation of the idea that the universe 
is on average homogeneous is a difficult one [6]. The 'solutions with isotropic expansion' 
referred to above are those where the evolution of the fiat metric on the torus is just a 
rescaling by a spatially constant conformal factor. If the expansion is not isotropic then 
certain parts of the geometry can be given arbitrarily and are not fixed by the field equa- 
tions. Thus, as we will see shortly, the field equations of the Newton-Cartan theory, unlike 
the Einstein equations, are not strong enough to determine a solution uniquely in terms 
of initial data. 

The results of this paper on the global dynamics of Newtonian cosmological models are 
linked to two questions of interest in cosmology. The first is the cosmic no-hair conjecture. 
This says that in a cosmological model with positive cosmological constant perturbations 
should be strongly damped. This conjecture refers originally to an infiationary period in 
general relativity. It is interesting to consider whether the analogous statement is true in 
Newtonian theory. In particular this may give insight into the mathematical mechanism 
responsible for the damping. Because of the quantities which can be given freely in the case 
where the expansion is not isotropic, it is clear that not all perturbations can be damped 
in the Newtonian case. This fact has been remarked upon by Barrow and Gotz [1] . (These 
authors did not use the Newton-Cartan theory.) Our results show that in the isotropic case 
small perturbations are strongly damped and at the same time set limits on how strong 
the damping can be. The second question is that of structure formation in cosmology. 
This has been studied in the context of Newtonian cosmology with a positive cosmological 
constant and dust as a matter model by Bildhauer, Buchert and Kasai[2]. They use the 
formation of singularities in a dust solution as a signal for structure formation in a more 
realistic model. Their results are consistent with those of the present paper and indicate 
that singularities are likely to occur in finite time for large initial perturbations. 

The paper is organized as follows. Section 2 describes the Newton-Cartan theory and sets 
up the basic equations needed in the sequel. In section 3 the results on global existence 
and asymptotic behaviour are proved. The examples which show the sharpness of these 
results are presented in section 4. 

2. Newtonian cosmology 

In this section we present some general information about spatially compact cosmological 
models in the Newtonian theory of gravity. More precisely, since Newtonian gravitation as 
it is normally presented is only applicable to isolated systems, we use a slight generalization 
of the standard Newtonian formalism, the Newton-Cartan theory. The original treatment 
of Newtonian cosmology, due to Heckmann and Schiicking [8] , is more elementary but does 
not bring out the underlying geometrical structures. The definition which will be used 
here is that which has been discussed by Ehlers[5] in the context of the Newtonian limit of 
general relativity. The basic objects of the theory are a covariant symmetric tensor gap (the 
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time metric), a contravariant tensor h°'l^ (the space metric), a torsion- free connection F^^ 
and a symmetric tensor T°'^ (the matter tensor) defined on a four-dimensional manifold 
(spacetime). They are supposed to satisfy the following conditions: 

1) At each point of spacetime there exists at least one vector with gaf3V°'V/^ > 0. 
Vectors with this property are called timelikc. 

2) If is timelike then the restriction of h"'^ to the space of covectors cja with (jUaV"' = 
is positive definite. 

3) gaph^'' = 0. 

4) gai3 and h°'^ are covariantly constant with respect to T^^. 

5) The curvature tensor R'^^s symmetry property h^'^Rp^^ = h'^^'RJaid- 

6) Rcp = 871 (g^a 9 PtT'''' - \g^TT'''' gad) - ^gaf3 where Rai3 is the Ricci tensor of T^^ and 
A is the cosmological constant. 

7) VaT«^ = 0. 

It is a known consequence of these conditions that the 3-dimensional distribution defined 
as the kernel of g^p is integrable. The integral manifolds are the surfaces of constant 
absolute time t. A spatially compact cosmological model in this theory can be defined 
as one where the hypersurfaces of constant time are compact. The tensor h°'f^ defines a 
Riemannian metric on each time slice and the field equation (condition 6) above) implies 
that this metric is flat. This means in particular that if a time slice is compact it must be 
isometric to a quotient of a flat torus by a flnite group of isometries[12]. 

By analogy with general relativity, a Newtonian spacetime M will be called time-orientable 
if there exists a global smooth timelike vector field V". It can be assumed without loss of 
generality that gapV^y^ — 1- Let S be one of the time slices and define a smooth mapping 
/ from an open subset U of x R to M by the condition that /(x, s) is the point of M 
obtained by starting at x and following an integral curve of V"' for a parameter distance 
s. The set U is supposed to be the maximal one for which such a mapping is defined i.e. 
each integral curve is followed to the end. The normalization condition ensures that the 
hypersurfaces t=const. coincide with the images of the hypersurfaces s=const. under /. 
This in turn implies that U is of the form S x I for some open interval /, since otherwise 
the time slices would fail to be compact. The mapping / is a local diffeomorphism. In 
particular f{U) is open. This construction can be started from any slice and so if / were 
not onto it would be possible to write M as a disjoint union of non-empty open sets. If 
spacetime is assumed to be connected this is impossible and so / is surjective. To get 
the injectivity of / it seems to be necessary to make another causality assumption. Call 
a subset of a Newtonian spacetime achronal if there exists no timelike curve starting and 
ending there. The following lemma can now be stated. 

Lemma 2.1 

Let M be a Newtonian spacetime with the following properties: 
(i) it is time orientable 
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(ii) each time slice is compact and achronal 



Then M is diffeomorphic to S x I, where S is a compact manifold admitting a flat metric 
and I is an open interval. 

As already remarked, such a manifold S is covered by a torus. Passing to this covering 
does not change the dynamics and so when investigating dynamical questions it may be 
assumed without loss of generality that 5" is a torus. This will be done in the following. 

Represent S in the form R*^/ ~, where the projection from to S is smooth. The tensor 
hab on M defines a tensor on R^ x / which will be denoted in the same way. Let habit) be 
its restriction to R^ x {t}. This is a complete flat metric on R^ and so must be isometric 
to the standard metric on R^[12]. This implies that given any orthonormal frame at the 
origin there exist coordinates on R^ x {t} such that the metric takes the form 5ab and the 
given frame takes the form d/dx"'. Introduce coordinates of this type on each slice. Since 
these can be constructed by following geodesies the are smooth functions on R^ x /. In 
these coordinates the components hab = ^ab- The identification which leads from R^ to S 
is of the form that x°' is identified with a;" + n*e^(t) for each triple of integers n*. What is 
not a priori clear is that e"(t) is smooth (or even continuous). A covering transformation 
of R"^ X / is smooth and must be of the form i— » + ?i*e"(t) for each fixed t. This 
is only possible if e"(t) is smooth. Doing a time dependent linear tranformation of the 
coordinates leads to a coordinate system where the identifications are time independent 
and the components hab only depend on t. 

The picture is now that the spacetime is the product of a torus with an interval and there 
are periodic coordinates such that the time and space metrics take the form 



Certain restrictions on follow from (2.1) and assumption 4) above. They are 



Names will now be assigned to the parts of the connection which do not vanish as a 
consequence of (2.2). 



fi'OO — I7 gOa — 
^Oa ^ ^ab ^ 




(2.1) 




(2.2) 




(2.3) 
(2.4) 
(2.5) 



Combining (2.2) and (2.4) gives 



(2.6) 
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This means in particular that &ab only depends on t and not on the spatial coordinates. 
Let p = T°°. Condition 6) gives the equations 



V„£;" = h^'dtOab - O'^'Oab - ^'''^ab + 47rp - A, (2.7) 

h'^^Vcnab = 0, (2.8) 

where indices have been raised with h"'^ while condition 5) implies the equations 

V[,Qab] = 0, (2.9) 

dtftab^VlbEa]. (2.10) 

Combining (2.8) and (2.9) shows that Qab is also independent of the spatial coordinates. 
Integrating (2.10) over the torus and applying Stokes' theorem gives 

dtQab = 0. (2.11) 
Substituting this back into (2.10) shows that the one-form Ea is closed. Thus 

Ea = Va- ^aU, (2.12) 

where r}a depends only on time and the function t/, which is determined only up to addition 
of a function of t by (2.12), is supposed to satisfy J U = 0. Consider now the effect of 
doing a coordinate transformation of the form x'" = + y", t' = t where the y"- depend 
only on t. This preserves the conditions (2.1) and changes Ea by habX^- Thus it can be 
used to eliminate 77°' without disturbing the assumptions made up to now. It follows that 
it can be assumed without loss of generality that 7]°' = 0. Before (2.7) is analysed ©06 will 
be split into trace and tracefree parts: 

©ab = Sob + l^hab, (2-13) 

where h'^^Eab = and ^ = h'^'^Qab- Then (2.7) becomes 

AhU ^-[9 + ^9^ + Eab^'''' - ^ab^"'' + 47rp - A] (2.14) 

where AhU = N^^VaS^bU. Let p{t) denote the average density at time t, i.e. p{t) = 
f pit, x)dx/ f dx. Then averaging (2.14) over the torus and using the fact that 6, Hab and 
Qab are spatially constant gives 

d+le^ + T.abT."^ - ^ab^"^ + 47rp - A = 0. (2.15) 

Putting this back in (2.14) shows that 

A^C/ = -47r(p-p). (2.16) 

Thus the Poisson equation has been recovered in a familiar form. 
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It will now be discussed how the solutions of the equations can usefully be parametrized. 
Taking the zeroth component of the equation VqjT'*^ = of condition 7) and averaging 
over the torus gives 

dtp + ep^Q. (2.17) 

To specify a solution we need to give the following objects: the matter variables (which in 
particular allow p to be calculated), 6 and Vtab on the initial hypersurface t = and Tiab 
everywhere. Then, using (2.6), (2.11), (2.15) and (2.17) the quantities hah, P, 9 and VLah 
are determined everywhere. All Christoffel symbols can be calculated from hah-, 9, 'Eab, 
ilab and VaU. The only one of these which is not already known is VaU. For reasonable 
matter fields the system of matter fields coupled to (2.16) should have a well posed initial 
value problem. In the case of a perfect fluid the equations look as follows: 

dtp + Va{pu'') + dp = (2.18) 
dtu"" + u^Vbu"" + p-'^h'^^Vw + l^w" - V'U + 2(S^ + ^Du^ = (2.19) 

assuming that the density p does not vanish anywhere. These equations are obtained by 
choosing the matter tensor 

T'^^ = pu'^u^ + ph'^^ (2.20) 

with vP = 1 and computing the condition 7) explicitly. That this procedure does in fact 
work for a fluid will be shown in the next section (Theorem 3.1). 

The simplest solutions of these equations are the Friedmann-like solutions which are ob- 
tained by taking p = p, u°' — Ea — 0, T^ab = and Q,ab = 0. In this case (2.6) reduces 
to dthah = ^Ohab, which implies that hab = cb'^{t)hab{0) for some function a{t). Suppose 
without loss of generality that a(0) = 1. This function satisfies the equation d = ^9a. 
Substituting this into (2.17) gives 

pa^ = p(0). (2.21) 
Combining (2.15), (2.17) and (2.21) results in the equation 

d = -^p(0)a^2 + iAa. (2.22) 

For solutions of (2.22) the quantity E = — ^p{0)a~^ — |Aa^ is constant and so the 
qualitative behaviour of its solution can easily be analysed using an effective potential. 
(The equation (2.22) also describes homogeneous and isotropic dust solutions in general 
relativity and in that case E = —k/2 where k is the curvature of the space sections.) This 
analysis shows that in the case A > 0, which is the one of greatest interest in the following, 
a tends to infinity as t — > oo precisely in the cases: 

(i) E > -|(47r)2/3AV3p(o)2/3 

(ii) E < -i(47r)2/3AV3p(o)2/3 and a(0) > (47rp(0)/A)i/3 

The constancy of E implies that for these solutions a/a ^ -^/A/B as t — > oo so that 
asymptotically the expansion becomes exponential. 
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3. The main theorems 



Now general solutions of the equations of Newtonian cosmology of the previous section 
with a perfect fluid as a matter model will be considered. The first question to be asked 
is whether the system of equations consisting of (2.16), (2.18) and (2.19) has a well posed 
initial value problem. There are various ways in which (2.18) and (2.19) can be written 
as a symmetric hyperbolic system for fixed U if an equation of state p — f{p) with /' > 
is given. In the following only polytropic equations of state will be considered i.e. p = 

71+1 

Kp^^ , where K and n are positive real numbers. It will be convenient to use the variable 
w = y/2n{2n + 2)Kp^/2n introduced by Makino [14]. Equations (2.18) and (2.19) can then 
be written as follows. 

dtW + u''daW+^wdaU'' + ^ew = (3.1) 

dtu"" + u^dbu" + ^wh^^dbw - h^^VbU + f^w" + 2(S% + n%)w'' = (3.2) 

The term in this system involving U is —V^U and this can be rewritten formally as 
47rV"A^^(p — p). This should be thought of as a lower order term, since it is formally 
of order -1. Accepting this for the moment, we see that contracting (3.2) with hac and 
taking the result together with (3.1) gives a system which is exactly like a symmetric 
hyperbolic system except that it contains one lower order term which is non-local. It is 
then natural to prove local in time existence for (2.16), (2.18) and (2.19) by following the 
steps in the standard local existence proof for symmetric hyperbolic systems as given in 
[13] for instance and checking that the non-local term does not cause any problems. This 
has been done in [3] for the case Tj^b — 0, Qab — 0. The general case is not significantly 
different. There results the following theorem: 

Theorem 3.1 

Let Tiab be a given C°° function oft. Let he an antisymmetric matrix, hP^^ a symmetric 
matrix, 9q a real number and po, Uq functions in H'^{T^) for some s > 3. Suppose that 
p is positive. Then there exists a T > and a unique solution of equations (2.6) (2.11), 
(2.15), (2.16), (2.18) and (2.19) on [0, T) such that 

(i) 0.ab, hab o-nd 9 are C°° functions oft which take on the values ^^b^ ^ab ^'^^ respec- 
tively for t=0 

(ii) p and u"" are C° with values in H'^{T^) and with values in H^~^{T^) and take on 
the values po and Uq respectively for t=0 

Moreover for this solution: 

(Hi) U is with values in H'^~^^{T^) and with values in H^(T^^) 

(iv) if p is uniformly bounded away from zero and the norms of p and u°' are bounded 

on the interval [0,T) then the solution can be expended to an interval [0,Ti) with Ti > T 

Next the global behaviour of solutions with almost homogeneous initial data will be dis- 
cussed. At this stage it will be assumed that /lab(O) = Sab^ "^ab = and O^b = 0. 
Some remarks will be made later on what happens when these quantities are allowed 
to be non-zero. Under these restrictions the spatial metric can be written in the form 
hab{t) = a'^{t)5ab- Then 6 = Sd/a. The equations (3.1) and (3.2) simplify and it is possible 
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to obtain a symmetric hyperbolic system from them in a shghtly different way from what 
was done previously. If v"" = au"' is taken as a new variable then the system is symmetric 
hyperbolic without contracting (3.2) with the metric. This has the advantage that the 
coefficient matrix of the time derivatives in the resulting system is the identity and this 
simplifies the derivation of energy estimates. Writing V for the pair {w — iv, v°'), the system 
is now of the form 

dtV + A^daV + BV + PV = (3.3) 

Here B = diag{3d/2na, d/a, a/a, a/a} and PV = (0,aPiV), where PiV is the gradient 
of the potential U of the density corresponding to w. All the matrices A'^ have a similar 
form. For example 

/ v'^/a w/2na 



A' = 



w/2na v^/a 

t'V«- 

V v^/a. 



In deriving estimates for these equations we will repeatedly use the following Moser in- 
equalities (see [11] for a derivation). 

P"(/^7)||2 < C(||/|U||i^^^7l|2 + PVIblblloo) 

WD^ifg) - fD^gh < C{\\Df\U\D^-'g\\2 + PVIblblU) 
IP"(i^(/))ll2 < C\dF/df\c.-4fro^'\\Df\U\D^fh 

Here the subscripts 2 and oo refer to the and L°° norms respectively, a is a multi-index 
and s = \a\. 

The basic energy estimate for (3.3) goes as follows: 
^11^112 = 2 J{V,dtV)dx 

= -2 J {V,A'd,V + BV + PV)dx (^-^^ 
< mA'Woo - 2min(|^, l)f )||K||i + 2||V^||2||P1/||2 
Now apply a spatial derivative D°' {a a multi-index) to (3.3) to get 

dtiD'^V) + A'diiD^'V) + BD^'V + D'^PV = -[D'^{A'diV) - A'D^'idiV)] (3.5) 
Applying the Moser estimates to the term on the right hand side gives 

\\D^{A\V)diV) - A\V)D^{diV)\\2 < C\dA/dV\cs' (1 + \\Vr-^)\\DV\\oo\\D'V\\2 (3.6) 
where s = \a\ and s' = max(l, s — 1). Thus 

f^WD^Vg < [C\dA/dV\cA^ + - 2min(A, 1) f] ^^^^^ 

+ 2\\D'V\\2\\D'{PV)\\2. 
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It remains to estimate PV. this will be done not in terms of w but in terms of p — p. 
Standard elliptic theory gives the estimate 



\\D^{PV)\\2<Ca\\p-p\\Hs-i (3.8) 

This inequality is only helpful if we have an estimate for p — p independent of (3.7). To 
obtain this, write the continuity equation in terms of a^{p — p), where k is some number 
< 3 to be fixed later. The result is 

dt{a\p - p)) + af^u'Viip - P) - '^^^a^p - p) = ^pW^u' (3.9) 

a 

Multiplying this by a^[p — p) gives an energy type estimate for p — p. The equation (3.9) 
can be differentiated in space to get higher energy estimates. The final result is 

j^\\a'D^-\p -p)\\l< {C\\DuU + ^^)\\a^D^-\p - p)\\l ^^^^^^ 



+ C[\\D{a\p - p))||oo + \\a\p - p)||oo + \\^p\\oo]\\a''D^-\p - p)\\2Wu\ 



2- 



The inequalities (3.7) and (3.10) are almost what is needed for the global existence theorem 
but still need to be altered slightly. First note that the expression C\dA/ dV\Qs' {'^ + 
11^11^^) can be thought of as a constant depending continuously on ||V^l|oc- Next, the 
dynamics of the scale factor a implies that there exists a constant Hq > 0, only depending 
on the initial data, such that a/a > Hq as long as a solution exists. It is also obviously 
possible to replace LP' norms of derivatives of order s of F and PV by norms of these 
functions themselves. The norm of PV can be eliminated in favour of p — p using (3.8). 
Finally the whole inequality can be divided by Hence 

^II^IIh^ < [C(||F|U)||W|U-niin(^,l)iyo]||F||H«+Cai-'^||a'=(p-p)||H.-i (3.11) 

In (3.10) the quantity Hq can be introduced in the same way. The inequality can be 
divided by ||a'^(p — p)||^s-i, giving 

j^\\a\p - p)U^-. < [C\\DV\\oo + {k- 3)i/o]||a'=(p - p)\\h^-. ^^^^^^ 

+ \\V\\Hs[\\Dia''ip - p))||oo + Wa'ip - p)||oo + llaVlloo]. 

The inequalities (3.11) and (3.12) are the main tools needed to prove the main theorem. 
Theorem 3.2 

Let initial data be given for (3.3) as in Theorem 3.1 with s > 4. Lf ao, 9q, po o.nd A > 
are fixed then there exists a positive number e such that if 

IIpo - Po||if« + WvqWh" < e (3.13) 
11 



the corresponding solution (whose local existence is guaranteed by Theorem 3.1) exists 
globally in time. Moreover the quantity \\v\\h'^ converges exponentially to zero and {p — p)/p 
tends to a limit in H^~^ as t ^ oo. 



Proof 

Choose some k lying strictly between 1 and 3 and define X — \\a'^{p — /5)\\hs-i + \\V\\h^- 
Let h = min((3 - k)Ho, 3/2n, 1). Then it follows from (3.11) and (3.12) that 



Here C{X) can be chosen to depend continuously on X (and hence to be bounded on 
bounded subsets) and / = max(l — k,k — 4). The assumption s > 4 is necessary so that 
the expression \\D{a'' {p — p))\\^ occurring in (3.12) can be estimated by \\a^{p — p)\\hs-i. 
The essential idea of the proof is now as follows. If X is sufficiently small then the first 
term on the right hand side of (3.14) will be negative and so will tend to cause X to 
become smaller. However more information is needed in order to conclude that X actually 
gets smaller since this term must be negative enough to outweigh the effect of the second 
term. This will be the case if t is large enough (and hence is small enough). Thus it is 
appropriate to do a bootstrap argument on a sufficiently long time interval. 

Let di be a positive number which is small enough so that X < di implies that C{X)X < 
h/3. If for a given solution of (3.3) on an interval [0, T) the quantity X is bounded then 
the solution can be extended to a larger time interval. Thus the inequality (3.14) shows 
that given any Ti > there is some ^2 > such that if X is less than 82 for t = the 
solution exists up to time Ti and X < Si on the whole interval [0, Ti] . Now choose Ti so 
that Co'(t) < h/3 for all t>Ti. (This is possible since a — cxo as t — > cxo.) Choose e > 
so that the inequality (3.13) implies that the initial value of X is smaller than 62- Now 
consider an initial datum for which (3.13) is satisfied. Let T* be the largest time (possibly 
infinite) for which the corresponding solution exists on [0,T*) and X < 5i for t < T*. 
Clearly T* > Ti. On the interval [Ti,T*) (3.14) implies that dX/dt < -^hX. It follows 
in particular that X is strictly decreasing there. Hence X(t) < X{Ti) on [Ti,T*). If T* 
were finite we could conclude firstly that the solution extends to a longer time interval 
(since X is bounded) and secondly that the inequality X{t) < 5i holds on a longer time 
interval. The second point follows from the fact that X{t) < X{Ti) < 61 on [Ti,T*) and 
continuity. This contradicts the definition of T*. The conclusion is that T* = 00. So 
the solution exists globally and dX/dt < —\hX for ever. A simple comparison argument 
then shows that X{t) < X(Ti) exp(— -i/it) for all t > Ti. Going back to the definitions 
shows that H^Hhs = is 0{exp{ — ^ht)) as t — > 00. Since a itself tends to infinity 

it follows that IIvHh" and HmUhs converge to zero exponentially as t — > 00. To get the 
remaining statement of the theorem, consider the inequality (3.12) in the case k = 3. 
Writing Y = \\{p — p)/ p)\\hs-i this leads to an estimate of the form 



dX 



< [C{X)X - h]X + Ca^X. 



(3.14) 



dY_ 
Itt 



< CX(Y + 1) 



(3.15) 
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Applying Gronwall's inequality and using the fact that X{t)dt < oo shows that Y is 
bounded. Using equation (3.9) in the case k = 3 gives the inequality 

p p Jtl 

for any two times ti and ^2 with ti < t2- This implies that there exists a function / such 
that {p — p)/p ^ f in H^~'^ as i ^ oo. ■ 

Theorems 3.1 and 3.2 concern fluids with polytropic equations of state. Very similar results 
can be proved for dust using the same methods. The differences will now be indicated 
briefly. Firstly, no Makino variable is necessary in the dust case; the density can be used 
directly. Secondly, the equations (2.18) and (2.19) are only coupled via the potential U 
when the pressure vanishes. Hence instead of (2.18) and (2.19) together forming a system 
which can be written in symmetric hyperbolic form, the equation (2.19) is symmetric 
hyperbolic by itself. Thus V can be taken to be just v"' in this case. Apart from this the 
argument works just as for a polytropic fluid. 

If is arbitrary then the analogue of the proof of Theorem 3.2 breaks down in general. 
However it is possible to generalize somewhat. Firstly it may be observed that flab makes 
no contribution in the energy estimates. Hence the analogue of Theorem 3.2 remains 
true if an arbitrary flab is allowed. In the case of non-zero T,ab it turns out that the 
condition needed to make the argument go through is that Qab should be positive definite. 
To ensure this it suflBces to ensure a sufficiently small uniform bound on S^. In general 
relativity it turns out that the shear, which is analogous to S^b, decays exponentially 
in homogeneous cosmological models with positive cosmological constant [18]. Thus the 
results found here can be regarded as indirect positive evidence for the validity of the 
cosmic no-hair conjecture in general relativity. 



4. Sharpness of the results 

In this section some examples will be presented which demonstrate the sharpness of the 
results obtained in the previous section. They all concern dust solutions and are proved 
by studying the paths of dust particles i.e. the integral curves of u"'. The solutions studied 
are such that p and only depend on t and one spatial coordinate, which will be denoted 
by X. Moreover it is asumed that the velocity is of the form ud/dx for some function u. 
Under these assumptions the equations reduce to 

Pt + {pu)^ + ep = o^ 

ut + uux + \eu - a~'^Ux = > (4.1) 
t/^rx = -47ra^(p-p) J 

The equation of motion of a dust particle is dx/dt = u{x{t),t). Differentiating this with 
respect to t and using (4.1) gives: 

d'^x/dt^ = i-pu + a-^Ua;)ix). (4.2) 
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Now consider two integral curves xi{t) and X2{t). 



^ix2 - Xi) ^ ~0^{X2- Xi) + I a '^U:e:rdx 

at 6 at Jxxit) 



3 dt 

Since pa^ is independent of time it follows that 



it) 

2 d /-^aCt) 
— 9—{x2 — xi) — Atx I {p — p)dx. 



(4.3) 



Xl{t) 



d n^'^ _ 3^ 3. , fdX2 dxi\ ^, 



On the other hand (4.1) and the relation 9 = 3d/a give 



r 
d/dt / 



pa-^dx = 0, (4.5) 

it) 



which just expresses the conservation of mass. Thus 

d'^/df{x2 - xi) = {Aa-^ - 2a/a)d/dt{x2 - xi) - Ba'^, (4.6) 

for some positive constants A and B. If d/dt{x2 — Xi) is initially non-positive then it 
will remain so. To see this note that if it is originally strictly negative it will remain 
so at least for a short time by continuity. If it is originally zero then (4.6) shows that 
it will immediately become negative and the same conclusion holds. Now suppose that 
d/dt{x2 — xi) becomes zero at some later time ti. Then the second derivative of X2 — xi 
at ti is non- negative and this contradicts (4.6). Thus d/dt{x2 — xi) remains negative as 
long as the solution exists. 

It will now be shown that for certain solutions of (4.1) the quantity (p — p)/p cannot 
converge uniformly to zero as t — > cxo. This shows that the fall-off property of the density 
perturbation proved in the last section is optimal. If the initial datum for p is not constant 
it is possible to choose starting values for xi and X2 so that 

/■ 0:2(0) 

/ p-pdx>0 (4.7) 

Now choose the initial velocity to be such that u{x2) —u{xi) < 0. Then u{x2{t)) — u{xi{t)) 
remains non-positive and so X2{t) — xi{t) < X2{0) — xi{0). Let 

M{t) = {a{t)f / p{t)dx = 2TT{a{t)fp{t) (4.8) 
Jo 
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Then M{t) is constant. Now 



/ (p - p)dx/p >C{ (p - p)dx)a- 

<-X2{t) 



rX2(t) 

= C[{ / pdx)a^ - ^{X2it) - xi{t))M{t)] (4.9) 
> C[i I ^ pdx)a^ - ^(a;2(0) - a;i(0))M(0)] 

Jxx(t) 

Hence J^^l^^^ {p — p)dx/p is bounded below by a positive constant and this is enough to 
show that {p — p)/p does not converge uniformly to zero. 

The argument just given applies for any value of the cosmo logical constant. It will now be 
shown that the analogues of the results of section 3 fail when the cosmological constant 
is zero. First the solution of (2.22) for A = and = will be computed. In that case 
a = ^87rp(0)/3a~-'-/^ and hence, after a translation in time, the solution is of the form 
a = Kt^/^. In particular the solution expands for ever. Now choose data once again for 
the fluid and for two dust particles so that dx2/dt — dxi/dt < as long as the solution 
exists. Then (4.6) implies that 

d^/dt^{x2 - xi) < -U~^d/dt(x2 - xi) - BK-H'^ (4.10) 



Consider now the corresponding differential equation 

d^y/df + It-^dy/dt = -BK-H'^ (4.11) 



Its general solution is 

y = -3BK-^logt-3Cit-^/^ +C2 (4.12) 

The qualitative behaviour of {xi —X2){t) can be studied by comparing it with the solution 
y{t) of (4.11) with the same initial data. Equation (4.11) is a first order ODE for dy/dt and 
so standard comparison arguments for solutions of first order ODE's give -^{x2—xi) < 
This inequality can then be integrated again to show that X2—X1 < y. Now (4.12) implies 
that J/ ^ in finite time if the solution of (4.1) exists that long. The trajectories of two 
dust particles in a regular solution of (4.1) cannot cross and so we have found data for (4.1) 
with the property that the corresponding solutions break down in finite time. Moreover 
these data can be arbitrarily close to homogeneity. 
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